Broad-band magnetic brain activity during rhythmic tapping tasks 
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1 Introduction 

Humans are able to generate and maintain rhythmic 
activity (as in running or beating a drum) and can also 
synchronize an internally generated rhythm with an 
external one. Brain functions involved in this ability 
include at least: internal rhythm generation, predic¬ 
tion of anticipated rhythmic stimuli, movement plan¬ 
ning and timing, and coordination of the brain areas 
responsible for these functions. Several magnetoen- 
cephalographic (MEG) studies have implicated var¬ 
ious brain areas in the ability to perform rhythmic 
tasks, including right frontal cortex [1], primary and 
secondary auditory cortices [2], and right frontal, pre¬ 
motor, and secondary auditory areas [3]. Functional 
magnetic resonance imaging (fMRI), positron emis¬ 
sion tomography (PET), and lesion studies [4, 5, 6] 
have also demonstrated the involvement of these ar¬ 
eas during timing and synchronization tasks. 

In this paper, we investigate the magnetic brain ac¬ 
tivity of subjects performing two rhythmic tapping 
and synchronization tasks. We use an adaptive linear- 
filtering technique known as differential synthetic- 
aperture magnetometry (dSAM) to locate sources of 
activity that have high signal-to-noise power ratios 
in various frequency bands. We then project the 
whole brain activity onto the individual sources with 
SAM spatial filters and examine the projections for 
evidence of phase-synchrony between the different 
sources. By isolating the oscillatory electromagnet¬ 
ic signals produced by the brain during performance, 
we can investigate the hypothesis that coordination 
between brain areas may be accomplished via phase¬ 
locking of the oscillatory activity. 

2 Methods 

2.1 Experiments 

Exp. 1: Missing Stimulus. Two subjects, one male and 
one female, were instructed to synchronize right (pre¬ 
ferred) index-finger flexions with an auditory metro¬ 
nome that had a fixed period of 800 msec (1.25 Hz). 
The metronome pulses were 80 msec 1 kHz tones. 


The 800 msec period has been shown to be very sta¬ 
ble in such synchronization tasks; even if a stimu¬ 
lus is omitted, subjects will continue tapping syn¬ 
chronously with the metronome [7]. In order to do 
this, subjects must predict the time of an upcoming 
metronome beat and plan a movement to occur at the 
same time. Brain structures that are responsible for 
the internal representation of the metronome rhythm 
and for motor planning must therefore be active dur¬ 
ing the task. We randomly omitted 12% of the stimuli 
in order to probe these brain structures. Subjects per¬ 
formed 21 blocks of 40 cycles each, yielding approx¬ 
imately 100 missing-stimulus responses. 

Exp. 2: Constant vs. Random. In this experiment, the 
metronome produced alternating 15-cycle trials of ei¬ 
ther constant 800 msec cycles or random 400-1200 
msec cycles with a mean of 800 msec. Three sub¬ 
jects (who did not participate in Exp. 1), two male 
and one female, were instructed to synchronize right 
index-finger flexions with the metronome during the 
constant period trials and to continue at the same tap¬ 
ping frequency during the random period trials. There 
were 16 trials (8 of each condition); the first three cy¬ 
cles of each trial were discarded to avoid transient ef¬ 
fects. When the metronome period was random, sub¬ 
jects were required to use their internal memory of the 
rhythm to produce constant period tapping. By com¬ 
paring the differences in brain activity between the 
two conditions, we can locate brain structures which 
are active during internal rhythm generation and self- 
paced movements, as well as areas related to exter¬ 
nal attention and synchronization. Auditory-only tri¬ 
als were also recorded in both experiments. 

Magnetic brain activity was recorded using a whole- 
head MEG system consisting of 148 radial magne¬ 
tometers (BTi Magnes). The data were sampled at a 
rate of 678 Hz using a high-pass filter with a low- 
frequency cutoff of 0.1 Hz. The MEG system was 
located in a magnetically shielded room (Siemens) 
placed on an active vibration-damping floor. MEG 
data were digitally filtered off-line in several frequen¬ 
cy bands including 15-25 Hz, 18-22 Hz, and 38-42 
Hz. Narrow bands are required for the phase analysis. 




2.2 Differential SAM 


be written as 


SAM filters are optimized to produce the highest 
possible signal-to-noise ratio for a given oscillatory 
source in the brain [8]. 



£ = h T X 
S 2 = h T Ch 

where C is the covariance matrix XX T of the band¬ 
pass filtered raw MEG data X (sensors x time), g 
is a column vector containing the theoretical magnet¬ 
ic field generated by a given fixed source, h is the 
SAM filter, £ is the resulting estimate of source ac¬ 
tivity, and S 2 is the source power estimate. The co- 
variance matrix was regularized by adding 5 times its 
smallest eigenvalue to each element of the diagonal. 
Differential SAM (dSAM) compares the level of 
source activity between active (a) and control (c) 
conditions, by computing the covariances and source 
power estimates separately and using the statistic p = 
s/Sj/S^. For large N, p is distributed approximately 
normally about a mean of 1. Sources with p values 
in the tails of this distribution are said to be differen¬ 
tially activated by the different conditions. We used a 
cutoff of p < .01 to determine active sources. 

2.3 Timeseries Analysis: Phase Locking 

The Hilbert transform x ^ z is defined as follows: 

2x2 if / > 0 
0 if / < 0 

z ■=$■ 2 

where =>■ represents the Fourier transform. 

If the Hilbert transform is applied to a bandpass- 
filtered real timeseries x, the result is a complex time- 
series 2 such that the real part is the original time- 
series x and the imaginary part is x phase-shifted by 
90°. The envelope of the original timeseries may 
therefore be calculated as \z\. 

Given the Hilbert transform z of a timeseries x, the 
complex phase 



is an estimate of the instantaneous phase of x at any 
given time. Given the complex phases from two dif¬ 
ferent timeseries x and y, their relative phase 6 may 



The complex relative phase 0 may be averaged over 
cycles to give the phase-locking value (PLV), a mea¬ 
sure of the extent to which the two timeseries are 
phase-locked [9]. 

«-v=$!|£«<o| 

The PLV for a given point in time is the length of 
the average relative phase vector in the complex plane 
(all the 9(i) have length 1). When two timeseries are 
phase-locked, the 6(i) will be approximately equal, 
and the PLV will approach 1. When the 9(i) are scat¬ 
tered randomly around the unit circle, the PLV will 
approach 0. 

In order to determine whether the PLV is statistically 
significant or not, we use a bootstrapping technique 
in which surrogate data are generated by shuffling the 
cycles of one of the timeseries [9]. If the two time- 
series are phase-locked within each cycle, but have 
varying absolute phases across cycles, then shuffling 
the cycles will destroy the phase-locking relationship 
and the PLV will be small. The phase-locking statis¬ 
tic (PLS) is defined as the standard deviation of the 
shuffled PLV values. The Hq -hypothesis that the two 
timeseries are not phase-locked may be rejected if the 
original (unshuffled) PLV is larger than some multiple 
of the PLS, typically 3. 

2.4 MRI-based Cortical Source Modeling 

SAM depends on an accurate current source model 
in order to calculate theoretical magnetic fields. We 
used the FreeSurfer package [10] to prepare white- 
matter surface models from structural MRI scans of 
the subject’s brains. 50,000 dipole sources were scat¬ 
tered randomly over the white-matter surfaces and 
oriented perpendicularly to them. Magnetic field- 
s were calculated using the homogeneous spherical 
conductor formula [11]. 

3 Results 

3.1 Behavior 

In Exp. 1, both subjects were able to produce 1-1 tap¬ 
ping. Mean Af ± SD (msec): Si, -186 ± 91; S 2 , 
— 194 Hf$4. In Exp. 2, one subject (Si) did not syn¬ 
chronize with the metronome at any time. The oth¬ 
er two were able to synchronize during the constant 
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Figure 1: Active sources (Si, Exp. 1) for the auditory - 
only vs. tapping contrasts in two frequency bands. All 
sources are significant at the p < .01 level, and in¬ 
clude left and right central and precentral sulci. 



Figure 2: Left: sources that were active during tap¬ 
ping, including left central and intraparietal sulci. 
Right: sources that were active for the missing stimu¬ 
lus, including left and right inferior frontal cortex. 

condition and reproduce the rhythm during the ran¬ 
dom condition. Mean IRI ± SD (msec): S 2 , 804 ± 45; 
S 3 , 806 ± 67 (constant); S 2 , 765 ± 73; S 3 , 814 ± 82 
(random). 

3.2 Differential SAM 

Fig. 1 shows primarily sensorimotor sources that were 
more active in the auditory-only condition than in 
the tapping condition, probably due to desynchroniza¬ 
tion during tapping. In Fig. 2 the missing stimulus 
activated BA44/45, while parietal and motor areas 
were desynchronized. Fig. 3 shows activation of right 
frontal cortex by the constant metronome and the left 
frontal cortex during internal rhythm generation (tap¬ 
ping) in the random condition. 

3.3 Phase Synchrony 

Both subjects in Exp. 1 had active sources in right su¬ 
perior frontal and left central cortices. The PLV val¬ 
ues between these sources are shown in Fig. 4. Be- 


Figure 3: Active sources for S 2 in Exp. 2. Left: con¬ 
stant vs. random in the auditory-only condition: right 
precentral and inferior frontal areas. Right: tapping 
vs. auditory-only during random trials: left inferior 
frontal and precentral areas. 

cause there were fewer cycles in the missing stimulus 
condition, the PLS values are much higher than for 
the normal condition, so it is unclear whether there 
was enhanced desynchronization prior to the move¬ 
ment or enhanced synchronization at —250 msec. 
Fig. 5 shows a constant vs. random contrast. Note the 
enhanced phase-locking in the random condition prior 
to movement. Also, the metronome had a desynchro¬ 
nizing effect (visible only in the constant condition 
since tapping was synchronized). 

4 Discussion 

We used rhythmic tapping tasks to establish patterns 
of coordination within the brain. By perturbing these 
patterns, we created contrasts that could be analyzed 
with differential SAM. We located sources in the right 
and left inferior frontal cortices that were activated by 
missing stimuli and internal rhythm generation, and 
which were phase-coupled with left precentral and 
central areas. 

Central and precentral areas were active bilaterally 
during unimanual tapping and exhibited temporally 
relevant variations in frequency and phase-synchrony 
with other regions, suggesting that ipsilateral central 
cortex may help mediate sensorimotor coordination. 
Additional active regions included occipital cortex 
(but not VI) and bilateral parietal regions (BA40/7, 
BA39), which also showed phase-locking with motor 
areas and may contribute to internal timing. 

Subject variability in patterns of inter-area phase¬ 
coupling may reflect differences in network organi¬ 
zation related to task strategies or performance skill. 
Frontal regions especially were variable from subject 
to subject which may indicate differences in internal 
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Figure 4: PLV values between right superior frontal 
and left central cortices in Exp. 1 (both subjects). Val¬ 
ues above 3 PLS are significant. Vertical dotted lines 
show flexion onset. 40 Hz band. 

rhythm generation strategies. 

Both event-related synchronization and desynchro¬ 
nization were observed, depending on modality and 
frequency. Some areas displayed almost continuous 
significant phase locking with event-related modu¬ 
lation. Correlated oscillatory activity may therefore 
signal the binding of populations of neurons while 
changes in the observed phase-couplings would indi¬ 
cate the flow of information between brain areas. 
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